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ABSTRACT 

All evolved stars with masses M < 2M Q undergo an initiating off-center 
helium core flash in their M c R3 O.48M He core as they ascend the red giant 
branch (RGB). This off-center flash is the first of a few successive helium shell 
subflashes that remove the core electron degeneracy over 2 Myrs, converting the 
object into a He burning star. Though characterized by Thomas over 40 years 
ago, this core flash phase has yet to be observationally probed. Using the Modules 
for Experiments in Stellar Astrophysics (MESA) code, we show that red giant 
asteroseismology enabled by space-based photometry (i.e. Kepler and CoRoT) 
can probe these stars during the flash. The rapid (< 10 5 yr) contraction of the 
red giant envelope after the initiating flash dramatically improves the coupling 
of the p-modes to the core g-modes, making the detection of £ = 1 mixed modes 
possible for these 2 Myrs. This duration implies that 1 in 35 stars near the red 
clump in the HR diagram will be in their core flash phase. During this time, the 
star has a g-mode period spacing of AP g « 70 — 100 s, lower than the AP g « 250 s 
of He burning stars in the red clump, but higher than the RGB stars at the same 
luminosity. This places them in an underpopulated part of the large frequency 
spacing (Au) vs. AP g diagram that should ease their identification amongst the 
thousands of observed red giants. 

Subject headings: stars: interiors — stars: oscillations — stars: late-type 



1. Introduction 



Sensitive space-based photometry from CoRoT (IBaglin et al. I l2009h and Kepler (IBorucki et al. 



20091 ) has enabled detection and characterizatio n of radial and non-radial acoustic (i.e. p 



mode) oscillations in thousands of red giant stars f|De Ridder et al.ll2009l : iBedding et al.ll2010 



Huber et al. 2010: Mosser et al. 2010; Hckkcr et al 



2011: 
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consistent with the scaling of the acoustic cutoff frequency (IBrown et al.lll99lf ). These modes 



generate photometric va riability ranging from 5 to 1000 pa rts per million (IMosser et al. 



2011bl : lHuber et al.lboilh and have lifetimes of at least 15 days ( jHuber et aliboiOt 



Baudin et al. 



2011). Their excitation and damping (and thereby their resulting amplitudes; iKjeldsen fc Bedding 
(199 5])) is related to the pres ence of vigorous convection within the star (see review by 
Christensen-Dalsgaardl (j2011al )). The majority of these oscillations are p-modes with nearly 
evenly spaced frequencies at 



Av Pa 4.27/iHz 
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a relation used with equation (TjQ) to determine the ste llar mass, M , and radius R from 
the measured v max , effective temperature T eff and Av flHekker et al] 1201 ll). The implied 
relation between v 



Hekker et al. 



2009 



2011b: Huber et al- 



and Av has also been confirmed observationally (ISte 



201ll ). 



Bedding et al. 2010; Huber et al. 2010; Hekker et al. 2011 



lo et al. 



2009; 



Mosser et al. 



Space-based observations ([Bedding et al.l |2010| ; iBeck et al] 1201 It IMosser et al] l2011a[ ) 
also enabled the detection of the angular degree £ = 1 mixed modes, which have p-m ode char 
acter i stics in the r e d giant envelope, but g - mode characteristic s in the helium core ( Scuflaire 



1974 



Osaki 



1975 



2004 : iDupret et al. 



2009 



Aizenman et al] 119771; Dziembowski et al] l2001t I Christ ensen- D alsgaar d 



Montalban et al. 



20101 ). Modes nearly evenl y spaced in period, at 



AP obs , around the I = 1 p-modes were identified by IBeck et~aL ( 2011 ) as a characteristic of 
the interior core g-modes allowing iBedding et al] ( 1201 ll ) to distinguish first ascent red giant 
branch (RGB) stars (i.e. those with degenerate helium cores) from red clump stars (i.e. those 
with non-degenerate He burning cores). This separation in the AP Q b s — Av diagram was 
also seen by CoRot (IMosser et al.ll2011aj ). and is a powerful new tool for stellar population 
studies. 

We show here that this new capability to probe the deep interior of a red giant should 
allow for the identification of those M < 2M stars undergoing the helium core flash. Known 
for more than 40 years as the defining e yent that ends the ascent of low mass stars up the 
RGB (and defines the tip of the RGB; ISalaris et al] (120021 )). this thermally unst able and 
off-center helium burning leads to a tf ~ 2 M yr phase of successive He subflashes ( iThomas 
19671 ; llben fc Renzinilll984t iMocak et al.ll2008h that remove electron degeneracy and convert 
the He co re to a stably bu r ning n on-degenerate object. However, there has been debate (see 
review by llben fc Renzinil (119841 ) ) as to whether the initiating flash can become dynamical 
in some way, or rather remains hydrostatic. 



Using the MESA code ( IPaxton et al.ll201ll ). we start in §2 by describing the changes in the 
red giant envelope and helium core during the core flash. We work in the Wentzel, Kramers, 
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Brillouin (WKB) approximation in §3, summarizing the p-mode and g-mode properties dur- 
ing the helium core flash. We close in §4 by discussing how stars in this phase of evolution 
can be differentiated from the more numerous populations of RGB and clump stars. 



2. Stellar Evolution during the Core Flash Phase 



Evolved stars with M < 2.OM undergo a He core flash that ends their RGB evo- 
lution, and leads them to the red clump where they undergo core He burning for « 70 
Myrs. The flash starts as an off-center thermally unstable ignition of helium burning at a 
mass coordinate ab out half-way through th e M n 0.48M^ He core ( Dominguez et al.ljl999j : 



S alar is et al. 



2002 



Paxton et al. 



Serenelli & Weissl bood : ISerenelli & Fukueital l2005[ iMocak et all bood ; 



20111 ). nearly independent of metallicity and M. 



Using MESA version 3709, we have calculated this evolution for a range of M, all with 
X = 0.7 and Z = 0.02, treating convection with the Schwarzshild criterion (i.e. no semi- 
convection, thermohaline mixing or convective overshoot) and a mixing length parameter 
of Qmlt = 2.0. We had no mass loss, diffusion, or rotation and we used the MESA "basic" 
(IPaxton et al.l 12011) nuc lear network with rates from NACRE. The opacity was from OPAL 



( Iglesias fc Rogers 



19961 ). with th e low temperature opacit ies taken from Ferguson et al. 



( 2005 ) with metal ratios given by iGrevesse fc Sauva l (1998) , as describ ed in IPaxton et al. 



(120111) . We used the OPAL eq uation-of-state ([Rogers & Navfonovi 120021 ). and HELM exten- 
sions (ITimmes fc SwestylbOOOl ) where needed. 



The evolutionary tracks in Figured] for a M = 1.0M Q and M = 1.8M star are for the 
time before, during, and after the core flash phase. The He core expansion from the initiating 
flash leads to an adiabatic temperature drop in the overlying H burning shell; quenching the 
burning. This loss of an energy source for the red giant envelope tri ggers a Kelvin- Helmholtz 
(KH) contraction from the tip of the RGB on a rapid timescale (IThomasI 119671 ). reaching 
L w 100 — 2OOL in 10 4 years (the open blue circles). 



Thomas! (119671 ) showed that the He core flash phase is governed by a series of shell 



subflashes that increase the entropy of the convective regions. This is seen in Figure [2j 
where five successive subflashes are seen to move inwards in mass coordinates. However, 
convective He burning is present for only about 10% of the time. The tj ~ 2 Myr timescale 
is set by the need for thermal d iffusion to act inwards between subflas hes, ending when the 



constant pressure (IPaxton et al 



thermal wave reac hes the core (Thomaslll967l ; ISerenelli &: Weissl 120051 ) . heating it at nearly 



20111 ) to the condition needed for the convective He burning 



core phase of the red clump. We end Figure [2] at that time. Even though the He burning 
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Fig. 1. — The evolution of a solar metallicity M = l.OM star and a M = 1.8M star before 
(black line), during and after (red line) the core flash phase. The filled circles (squares) 
show where the star resided on the RGB at 30 (60) Myrs prior to the initiating flash at 
L tip = 1O 3,4 L . The open circles show the location 0.01 Myr after the initiating flash, 
whereas the open blue square is 2 Myrs after the initiating flash, a time when the star 
reaches the red clump, where it still resides at 60 Myrs (green star). 
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lumino sities (blue dashed line s) get quite large (Lu e ~ 10 9 3 L Q in the initiating flash at 
t — 0), IShen fc Bildstenl (120091 ) have shown that the heating from He burning always occurs 
on a timescale much longer t han the loc al dynamical time, an outcome confirmed in multi- 
dimensions by iMocak et al.l ( 120081 . 120091 ) . The stellar response during the 2 Myr core flash 
phase is shown in Figure [2] and form the basis for our asteroseismological work. 



3. Mode Properties during the Core Flash Phase 



ied. 



The non-radial adiabatic mode structures for RGB and red clump st ars have been stud- 
with focus on the I = 1 mixed modes that reveal core properties (Dziembowski et al. 



2001 



2011a 



b; 



Jiang et al 



2011 



di Mauro et al 



Christensen-Dalsgaardl 2004 ; Dupret et al.l 2009 ; Montalban et al.lboiol ; Christensen-Dalsgaard 



20111 ). These mixed modes have a p-mode qual- 



ity in the outer parts of the star (where they are excited by convection, and typically have 
n n « 10 radial nodes), but penetrate into the stellar core as very high order {n„ > 100) K- 



modes. This justifies our initial exp loration in the WKB limit (jUnno et al.lll989l ; lAerts et al. 



20101 ; IChristensen-Dalsgaardll2011al ). where the local (at radius r) radial wavenumber, k r , is 



c 2 u 2 



(u, 2 - N 2 ) (c 2 - Sj) 



(3) 



where u = 2irv is the mode frequency, c s is the sound speed, N 2 is the Brunt- Vaisala 



frequency, and Sj 



+ l)/r 2 is the Lamb frequency. Figure [3] shows the propagation 



diagrams (radial profiles of N, Si and S2) for four phases of the M = l.OM model, all chosen 
when Au = 4^Hz. These resulted in nearly the same values of i/ max ~ 28/iHz (denoted by 
the horizontal line). From top to bottom, the panels are for t = —60 Myrs prior to the 
initiating flash (on the RGB), t = 1.4 Myrs (after the initiating flash), t = 1.6 Myrs (during 
a convective He burning subflash) and t = 2.0 Myrs when the star is in the convective He 
core burning phase of the red clump. 

Modes supported by the acoustic response of the envelope (i.e. p-modes) propagate 
in the outer parts of the star where u 2 > Sf and u 2 > N 2 . Since iV 2 is either zero or 
very small there, we set c 2 k 2 = uj 2 — S\. The eigenfrequencies are then found by setting 
J k r dr « n p 7r. In the extreme limit of u> 2 3> S 2 , this simplifies to n p 7r = u> J dr/c s , where the 
integral extends over the outer parts of the star where u 2 > Sf. Tradition is to write this as 
v = n p Au, where the "large spacing" Av is defined with the integral over the whole star, 



Av 



-1 



2/"-. 

JO c s 



(4) 



a reasonably accurate representation. 
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Time Since Initiating Helium Flash (10 6 years) 



Fig. 2. — Properties of a M = 1.0M Q star for the 2 Myr core flash phase. The blue shaded 
regions show where convection is present in mass coordinates, whereas the orange shaded 
regions are where nuclear burning yields > 1000 ergs gr _1 s _1 . The stellar luminosity is 
shown by the orange solid line, whereas the H (He) burning luminosities are shown by the 
green dashed (blue solid) lines. The stellar radius, R, is shown by the brown dot-dashed 
line, and the dotted brown line shows the He core radius. 
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Fig. 3. — Propagation diagrams for a M = l.OM star at four different stages of evolution 
(described in text), all when Av = 4/iHz. The Brunt- Vaisala frequency, N/2tt (solid lines), 
and Lamb frequencies Si/2n (dotted lines) and S2/27: (dashed lines) are shown. The thin 
horizontal line is at z/ max ~ 28/iHz. From top to bottom, these models have L/L Q = 
41.6, 50.6, 51.1&50.8, and R/R Q = 10.9, 11.2, 11.2&11.2. 
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Modes supported by the internal buoyancy (i.e. g-modes) propagate in the stellar inte- 
rior, where oo 2 < N 2 and oo 2 < S 2 . As evident from Figured! these inequalities are strong for 
frequencies ~ i/ max , yielding k 2 w £(£ + l)N 2 /r 2 uj 2 . The high radial order (n g 3> 1) g-modes 
then have n g 7r = (£(£ + l)) 1 ^ 2 / Ndkir/u, yielding v~ x = n g AP g , with AP g being the period 
spacing. For £ = 1, this is 

2 1 / 2 7T 2 

At z/ max , n g l/(z/ max AP 9 ) > 100, safely in the WKB limit. The integral J Adlnr extends 
over the region where u 2 < N 2 and u> 2 < S 2 . 

Now consider how AP g and Av evolve during the core flash phase. The KH contraction 
after the initiating flash triggers a rapid decrease to R ~ IIP©, and the increase of Av 
evident in the second panel of Figure HJ Thereafter, the value of Av simply responds to the 
mild radius changes during subsequent subflashes. The He core radius expansion and outer 
envelope contraction triggered by the initiating He flash leads to an outer envelope structure 
during the tf ~ 2 Myr that is nearly the same as that on the red clump. Indeed, the profiles 
of N 2 and S 2 in Figure |3]for the bottom three panels (all after the initiating flash) are nearly 
identical for r > O.O4P . 

The prime property that is changing over the 2 Myr is the He core, as it undergoes 
additional subflashes that lift the degeneracy of the deep interior. The £ = 1 g-mode period 
spacings for the models in Figure [3] are AP g = 61.2 s for the RGB, AP g = 73 s for the 
core flash phase model between subflashes (second panel), and AP g = 255 s for the red 
clump. During the convective He burning shell subflashes (less than 10% of the He core 
flash phase, so < 2 x 10 5 years), an intervening evanescent zone appears in the core (see 
third panel in Figure [3]). Though coupling through this zone is possible (see below), we 
calculated the period spacing assuming that only the outer g-mode cavity is relevant, giving 
AP g ps 250 — 270 s during each sub-flash. As evident in the third panel of Figure HI the 
existence of a degenerate core with N 2 > during the most of the core flash period keeps 
AP g < 100 s. The fully convective He burning core o f the red clump is what causes the 



increase to AP g = 255 s (1Christensen-Dalsgaardll2011at ) at the end of Figure HI 



The £ = 1 mixed modes are more likely to be detected when the coupling through the 
outer evanescent region (the region in the star where u 2 > N 2 and u 2 < Sf so that k 2 < 0) is 
strong. In the WKB limit, the ratio of the mode amplitudes (and also location of the mode 
energy) between the two turning points, the inner one at T\ (where u 2 < N 2 first occurs) and 
the outer one at r 2 (where u 2 > S 2 first occurs) is largely determined by / = k r dr, but 
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Time Since First Flash (Myrs) 



Fig. 4. — Seismic properties during the core flash, with vertical grey regions showing sub- 
flashes. From top to bottom, the panels show i/ max from equation ([I]), Av> from equation (j4]), 
AP g from equation ([5]), the Helium burning luminosity Ln e , and the evanescent integral /. 
The vertical dashed lines denote the times for the bottom three panels of Figure [3] 
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also depends on n and r 2 ( Aerts et al. 201ol )l*l The coupling of the I = 1 modes is stronger 
than that of £ = 2 simply b e cause the outer tu rning point is closer to the core for £ = 1 
( IChristensen-DalsgaardI 120041 : iDupret et al.ll2009l ). The values of I (calculated at u — ^ max ) 
for the models in Figure El are I = 1.26 for the RGB, and / « 0.3 for the other three models 
(which have nearly identical outer envelopes). The bottom panel of Figure H] shows that I 
changes little over the 2 Myr, giving us confidence that the seismic probe of the core will be 
as fruitful for the He core flash stars as the red clump. 

We have not fully investigated the mode structure during the subflashes. However, a 
few things can be said. The first is that an additional inner evanescent zone appears during 
subflashes (the region at 0.0115 < r/R Q < 0.0192 in the third panel of Figure |3]) splitting 
the g-mode cavity into two zones, each with their own distinct period spacing. This adds 
an additional evanescent integral Ij ~ [£(£+ l)] 1 / 2 cilogr « 0.72 that must be accounted for 
in the observable mode structure. Due to that extra penalty, we chose to display (in Figure 
|4| the period spacing during the subflashes as that of the outermost g-mode cavity. If the 
inner cavity proves to be adequately coupled, these extra modes would cause oscillations in 
the observed period spacing diagram. That full analysis awaits our future efforts. 



4. Detecting the Core Flash Phase 



Our re sults are best seen in r elation to the observational work performed by lBedding et al 



( 1201 ll ) and iMosser et al.l ( l2011al ). Their secure detection of mixed modes and measurement 
of their period spacing, AP obs , led to a clear distinction (see Figure |5]) of RGB stars (blue 
points) from the red clump stars (red points) of M < 2M Q and the extended red clump 
(yellow points). The evolution of AP g for our 1.0M Q (green lines) and a 1.8M (black lines) 
models are shown in Figure [5] as a dashed line for the RGB, a solid line during the core flash 
and a short-dashed line for the core He burning phase. 



Bedding et al.l (120111 ) discussed the challenge of secur e ly inferring AP g when only a fe w 



mixed modes are observed (see also IMosser et al.l (j2011al ); IChristensen-DalsgaardI ( 1201 lbl )). 
In particular, they found that the measured spacing, AP b s , would be less than AP g by 
a factor of 1.3 — 1.6. Such a need for scaling is evident along the RGB, and will be best 
resolved by longer duration Kepler data that reveals addi tional mixed modes . This should 
prove possible given their much longer expected lifetimes (IDupret et al.ll2009l ). 



Stars undergoing the core flash occupy a distinct region in Figure [5j The only other 



^dn the plane-parallel limit (e.g. T2 —T\ Crj), the mode amplitude ratio is oc exp(±/). 
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Fig. 5. — Evolution of l.OM (green) and 1.8M (black) red giants through the helium core 
flash in Av — AP g space. The RGB of a 2.4M star is shown by an orange line. The dashed 
lines are the RGB phase, whereas the solid line is the core flash (with the rare excursions 
during subflash es omitted) and short -dashed the He burning stage. The colored data points 
are AP Q b s from iBedding et al.l (l201lh . 
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models that traverse this region are M > 2.OM RGB stars, and potentially AGB stars, both 
reasonably rare instances that will need to be observationally distinguished. The simplest 
technique is to use the measured z/ max values, as z/ max ~ 38 Hz for the 2.2M and 2.4M Q 
RGB models, whereas z/ max < 32 Hz during the core flash for the l.OM model of Figure 
HJ All objects with 70 s < AP g < 200 s and 2/iHz < Au < 6/xHz are worthy of a serious 
analysis to see if they are explained as ~ 1M core flash stars rather than > 2M RGB 
stars. We highlight three such candidates in Figure [5] by slightly enlarging their blue points. 
Since AP t, s is always less than AP g , these points could only move upwards in this diagram, 
placing them securely in the region of interest. There are 193 red clump stars in this figure, 
so we exp ect about 193/35 ~ 5 flashing stars, consistent with the few "outliers" in this early 



data set. iMiglio et al.l ( 120091 ) showed that the expected metallicity of most of the CoRoT 
red giants is near solar, with a few as low as one-tenth solar. Our MESA runs for one-tenth 
solar occupy a comparable part of this diagram, though they go to slightly higher values of 
Az/ during the flash due to reaching a lower luminosity during KH contraction. 

There is much to be learned from the asteroseismic discovery and analysis of a star 
undergoing the core flash phase. It would be an immediate confirmation of the ~ 2 Myr 
duration of the event; eliminating the alternate possibilities of more dynamical outcomes. 
Study and analysis of the g-modes may reveal the specific stage of evolution for such a 
star in the 2 Myr phase. Any possible inferences regarding the rotation rate of the core 
during these 2 Myrs would be especially valuable, as the short timescales may not allow for 
perfectly rigid core r otatio n. Another path to detection of this rare phase was highlighted by 



Silva Aguirre et al.l ( 120081 ). who noted that RGB stars with substantial mass loss will cross 
the RR Lyrae instability strip on their way to the horizontal branch. Those stars would then 
exhibit a detectable pulsation period change due to their secular evolution. 

Additional calculations are certainly needed, both to fully explore lower metallicities, 
as well as to see how additional mixing mechanisms (e.g. thermohaline) may impact the 
qualitative evolution. Fully consistent pulsation calculations are needed, both to confirm 
the WKB estimates here, but also to provide damping times and to explore the interesting 
periods of subflashes (only relevant for f« 2 x 10 5 years) that may prove detectable if the 
sample size of red giants with measured AP g gets large enough. 



We thank T. Bedding, T. Brown, J. Christensen-Dalsgaard and V. Silva Aguirre for 
discussions. This work was supported by the NSF through grants PHY 05-51164 and AST 
11-09174 and by support for K.M and P. J.M. through the Worster Family Summer Fellowship 
fund. 
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